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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

CALCULATION OF THE EFFECTS OF STRUCTURAL FLEXIBILITY 
ON LATERAL CONTROL OF WINGS OF ARBITRARY 
PLAN FORM AND STIFFNESS 
By Franklin W. Diederich 

SUMMARY 


A method is presented for calculating the effectiveness and 
reversal of lateral- control devices on wings of arbitrary plan form 
and stiffness. Computing forms and an illustrative example are 
included. 

• The margin against aileron reversal is sho^m to he relatively low 
for swept wings at all speeds and for all configurations at supersonic 
speeds; the margin is relatively high at subsonic speeds. Effectiveness 
of conventional aileron configurations on sweptback wings at supersonic 
speeds is relatively low. 


INTRODUCTION 


Adequate lateral control constitutes one of the more important 
design requirements for airplanes. The ability of the airplane to 
enter a roll is determined by the control power and is measured by the 
maximum available rolling moment resulting from lateral- control deflection 
The degree of lateral maneuverability may be represented by the helix 
angle at the wing tips corresponding to the highest rate of roll. The 
lateral maneuverability depends both on the control power and the damping 
in roll . 

The control power and the damping in roll 'are affected by structural 
flexibility. Control deflection ordinarily gives rise to aerodynamic 
loads which tend to deform the wing structure in such a w^ as to reduce 
the loads on it and thus to reduce the control power. If the dynamic 
pressure of the air stream is sufficiently high, the amount of load 
which results from the structural deformation may be sufficient to 
nullify the effect of the control deflection. The speed and dynamic 
pressure corresponding to this condition are known as the lateral- 
control reversal speed or reversal dynamic pressure, since at a 
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slightly higher dynamic pressure a control deflection in a given 
direction would result in a rolling moment in a direction opposite 
to that of the moment on a similar rigid wing. 

Much of the early work on lateral- control reversal and loss of 
control due to structural deformations was done in Great Britain. The 
first published account of an investigation concerned with aileron 
reversal appears to he reference 1. Even at that early date, aerodynamic 
induction was taken into^account, hut an arbitrary wing-deflection mode 
was assumed. The lateral- control power was analyzed on the basis of 
the same assumptions in reference 2. The work done in Great Britain 
subsequent to the publication of these two papers has been concerned 
with more refined means of accounting for the actual stiffness 
distributions and for aerodynamic induction; at the same time a great 
deal of attention has been devoted to the simplification of the 
numerical work required to obtain results in practical cases. 

Work done on the problem of loss of lateral control due to wing 
flexibility in this country (references 3, k, and 5, for instance) 
now represents the same stage of development as British work in the 
field. Eeference. 3 presents convenient methods for determining the 
aileron-reversal speed and other related critical speeds of wings of 
arbitrary stiffness distribution; aerodynamic induction is taken into * 
account by means of ah over-all correction. Eeference 4 is concerned 
with the determination of the lateral maneuverability and control 
effectiveness. Aerodynamic induction la taken into account approxi- 
mately, and the method is applicable to wings of arbitrary stiffness 
distribution. The numerical work required for the analysis is 
fairly extensive, however. Eeference 5 shows a method for calculating 
the reversal speed by matrix iteration; the method is convenient and 
applicable to arbitrary stiffness distributions, but the integrating 
matrices are only approximate and aerodynamic induction is taken into 
accoimt only by means of an over-all correction, unless suitable influence 
coefficient matrices are used in conjunction with the method. 

.Although the foregoing methods and similar British work are 
generally satisfactory for calculating the lateral-control effective- 
ness and aileron-reversal speed of straight wings, they are inapplicable 
to swept wings. The present paper is concerned with an analysis of 
these problems for wings of arbitrary plan form, including swept plan 
forms, as well as arbitrary stiffness. The method is based on the- 
analysis of the loading of flexible wings presented in reference 6. 

Since suitable aerodynamic -influence coefficients are not yet available, 
aerodynamic induction is taken into account only as an over-all 
correction and a slight reduction of the load at the tip, as in 
reference 6. G?he method is formulated in such a manner, however, that 
aerodynamic-influence coefficients may easily be included as soon as 
they become available. 
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The numerical analysis required in any given practical case 
constitutes an extension of the calculations outlined in reference 6. 
Computing forms for the additional calculations required for an analysis 
of lateral- control effectiveness or reversal are presented in this 
paper. Their use is descrlhed in the section "Application of the Method." 
This section may he read without reference to the derivation of the 
method. As an example illustrating the method, the lateral-control 
effectiveness and reversal of the wing considered in reference 6 are 
a na lyzed in this paper. The reversal speeds of several wings derived 
from this wing hy shifting the elastic axis and rotating the wing are 
calculated to demonstrate some general effects of sweep on the aileron- 
reversal speed. 


SIMBOIS 


The symbols used in the analysis are those of reference 6 with the 
following additions; 


[A] 

[ae] 


CPg 


®2 



yi 


yo 

“6 

6 


e 


auxiliaiy aeroelastic matrix 
reversal matrix 

center of pressure of the load produced hy aileron deflection, 
fraction of chord from leading edge 

distance from the reference axis to the center of pressure 
of the load due to aileron deflection (positive rearward), 
fraction of chord 

matrix defined in equation (U) 

distributed torque, inch-pounds per inch 

lateral ordinate of inboard end of aileron, inches 

lateral ordinate of outboard end of aileron, inches 

angle of attack equivalent to unit aileron deflection 

aileron deflection measured in planes parallel to the 
direction of fli^t, radians 

moment- arm ratio 
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DERIVATION OF THE METHOD 
Assumptions 


The assumptions made in the following analysis are the same as 
those made in reference 6. In addition it is assumed that the angle 
between the aileron and the wing is constant eilong the span of the 
aileron. 


Air Loads 

The aerodynamic forces on a wing section with control neutral are 
given in reference 6. In keeping with the aerodynamic assumptions the 
loading due to the aileron deflection is considered to he the corres- 
pouding strip loading multiplied hy a reduction factor and rounded 
off at the tip. The effective section lift-curve slope appropriate to 
the aileron loading has approximately the same value as the effective 
slope appropriate to the linear- twist loading, as may he deduced hy 
comparing the different J values presented in references 2 and 3. 

The aerodynamic force on a wing section of unit width parallel to 
the direction of flight is then 


V 


= meq.c cos 


A (ots + 056) 


( 1 ) 


where is the aileron effectiveness factor (^C2/36)yr ^c^/^a) . 

(See fig. 1.) The angle of aileron deflection 6 is measured in 
planes parallel to the direction of fli^t. If it is desired to 
measure 6 in planes perpendicular to some reference line, for 
instance, the quarter-chord line, the value of Og is replaced hy a 
value . which is og multiplied hy the cosine of the sweep angle 

measured to the reference line. The moment of the running load about 
the elastic axis is 


t* = m^qeic^ cos A (ag - €056) (2) 

where e is the moment-arm ratio and 02© is the distance 

between the center of pressure of the load due to aileron deflection 
and the reference axis. (See fig. 1.) 
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The acc\imulated torque and moment referred to the elastic axis 
may then he written as • 


{t} = mgCiSAei^c^S 



( 3 ) 



Equations of Equilihidum ^ 

The equations of equilibrium may he set up and treated in the same 
manner as in reference 6. The result is. the relation 



where a and [aJ are defined in reference 6, and the auxiliary 
aeroelastic matrix A is defined hy 
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If structtiral influence coefficients of the type described in 
reference 6 are uaed instead of the stiffness curves, equation (5) 
tates the form 




(7) 


where a’ and [a’J are defined in reference 6 and 



0 

< 

1 


rO 1 

[a-] = 6i^Or[«B][K6] 

®lr\Cr/ 

[°] - N W 

c 

5 r 


( 8 ) 


Solution of the Equations 


The aerodynamic loading corresponding to a given aileron deflectioi 
may be obtained by writing equations (5) or (?) lo ilio following form: 


[l] - a[A] 


jttgj = - a[A]|ag6j 


(9) 


Once the right-^and side of the equation (9) is multiplied out, the 
twist distribution may be calculated by solving the simul- 

taneous equations of equation ( 9 )* The loading as well as the accumu- 
lated torques and moments may be obtained from eqiaatlons (l) to ( 4 ). 


In order to calculate the aileronr-reversal speed conveniently from 
equations (5) or (7), it is necessary to eliminate by expressing 


it in terms of 




}• 


The required relation results from the condition 


that at the reversal speed the, rolling moment yanlshes, so that 
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where L^2Ji first row of the [E2I niatrix. The solution of 

equation (lOa) may he shown to he (hy the reasoning of reference 5) 



where the matrix 



is defined hy 


(U) 




1 0 0 0 . . 

1 0 0 0 . . 

1 0 0 0 . . 

1 0 0 0 . . 


( 12 ) 


If eq.uation (U) is substituted in equation (5), the following relation ^ 
is obtained; 

.Tag]. = a [AE'|{a's} (I3) 

where the aileron- reversal matrix P^eI is defined hy 



The reversal dynamic pressure is calculated by iterating eq.uation (13) 
and substituting the critical value of a in the equation defining the 
parameter (eqviation (16) of reference S\. 
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APPLICATION OF THE METHOD 
Selection of the Parameters 


The section lift- curve 'slope and aerodynamic- center values are 
chosen for the Mach number of interest, as described in reference 6. 
At subsonic speeds the lift-curve slope is corrected for finite- 
span effects as described in said reference. Values of oq and cpg 

are best obtained from experimental section data at the appropriate 
Mach niamber. Theoretical thin-airfoil values of these parameters are 
presented in figure 1 for subsonic and supersonic speeds. 

The structural parameters are obtained in the manner described in 
reference 6. 


Calculation of Matrices 

Either a 6-point or a 10-point solution may be employed. Computing 
forms are provided for the 6-point solution j s imi lar forms may easily 
be set up for the 10-polnt solution. 

In order to tajce account of the location of the inboard and outboard 
eBctremlties of the aileron with the relatively few stations used in the 
analysis, equivalent 6 values have to be used. These values are 
given in figure 2. They are intended to give a rounded off 6 variation 
which has approximately the same area and the same moment about the 
root as the actual 6 variation. The equivalent 6 values of figure 2 
pertain to actual values of & equal to 1; they apply to ailerons which 
extend from yi/s^ combined to apply, to any 

aileron configuration. Several examples are listed below for the 
six-point method, the act'ual values of 5 being 1 and tlie equivalent 

y . 

values being read from figure 2(a) as O.716 for ~ = 0.55 and as 

yi 

0.293 for —=0.95: 

S-,.r 
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Case 

HI 

2 

3 

HH 

5 

(yi/3w) 

0.55 

0.95 


0 

0 

(yc/®w) 

1.00 

1.00 

•95 

1.00 

•55 

(y/®w) 


0 

0 

0 

0 

1 

1 

.20 

0 

0 

0 

1 

1 

.40 

0 

0 

G 

1 

1 

.60 

.716 

0 

.716 

1 

.284 

• 8o 

1 

0 

1 

1 

0 

• 90 

1 

.293 

.707 

1 

0 

(yi/sw) corresponds to the inboard 

extremity of the aileron or elevon 
(yo/®w) corresponds to the outboard 

extremity of the aileron or elevon 


G?he values for case 3 ar© obtained from those of cases 1 and 2, those 
for case 5 from the ones of cases 1 eind 4. 

The [a] matrix is calculated as described in reference 6. The 
calculation of the auxiliary aeroelastic matrix then proceeds as shovn 
in table Ij the munbering of the steps indicated in the upper left 
corner of each block are a continuation of the steps in table Vl(b) 
of reference 6. 

If it is desired to calcxilate the aileron -reversal dynamic pressure, 
the aileron -reversal matrix is calculated by means of the form of table II. 
The value of G is calciilatod by multiplying the first row of 
matrix C"®®® reference 6) by the column matrix In 

accordance with equation (l4), the -|ag6l values also occur as a diagonal 
matrix; in this form they are divided by the G value and entered in 
matrix [@] • The calculation then proceeds according to the 
instructions of table II . 
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Special cases arise when any or all of the ei or 09 values 
are zero. If only e^^ is zero, the e^^ .value at some other point 

may he used as a reference throughout the analysis and the parameter a 
redefined accordingly. The first column of the matrix Tfl^il is 


calculated in this case hy multiplying the first column of the 

®2 


matrix 


Ki 


hy the ratio 


value is zero, say the fourth along the span, ei is used as a 

j-r 


reference hut the fourth column of 
the fourth column of the matrix 


®lref. 
1© s 


Similarly, if some other 


[®] N 


is calcxilated hy multiplying 
©2 

hy — where 02,, v 

®1„ 


is the val,ue of eg at the fovirth station, at which e^^ is zero. 

If ej_ is zero along the entire span, some of the comtputihg 
instructions given in this paper, as well as the ones given in reference 6 


®2 /c 


®1 /c 


®2j.(cr) 

matrix is entered in the space provided for the 


V 


matrix. Some of the instructions of table VI (h) of reference 6 and 
table I of this paper are then modified as follows: 


Step ® 


Steps (8) to 


Step (1 


[Kl] 


®2rV°r. 


as indicated for the case e^ = 0 in reference 6. 

®2r®r ^ (El)j. 1 


(GJ)^ tan A 
All other Instructions are unaffected. 


If 62 is zero along the span, table I of this paper may he 
modified as follows;. 

Step (^) [a] = 

Step may he omitted; all steps in table II are unaffected. 
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Calculation of the Aileron-EereraaJ. Speed 

The [Aj^] matrix is iterated in table III (a) to calculate the critical 

veilue of the parameter a and hence the critical speed. The calculation 
ordinarily has to he performed at least twice, once for subsonic speeds 
and once for supersonic speeds. From' these critical values^ from the 
definition of the parameter a, and from the effective lift-curve slope 
the dynamic pressure required for aileron reversal may be calculated 

and plotted as a function of Mach number. If the actual dynamic pressure^ 
for the altitudes of Interest is also plotted on the same chart, the 
lowest intersection of the reversal with a true-dynamic-pressure line 
will give the reversal Mach number and dynamic pressure at the altitude 
of the true-dynamic-pressure line. 

The |Aj^^ matrices calculated for the special cases mentioned in 

the preceding section do not all yield the critical value of the 
parameter a. When the value of e^ is zero at the root, the critical 

value of the parameter a based on the reference value of e^^ will 
be obtained. If ej_ is zero at some other point along the span, or if 
02 is zero along the entire span, critical values of the parameter a 
will be obtained. In the case where e^^ is zero along the entire 
span, iteration of the Aj^ matrix calculated by following the instructions 

of the preceding section will yield the value of the parameter d at 
divergence . 

In some' of these special cases, and possibly in other cases as well, 
it may be found that the iteration procedure does not converge. In 
those cases th^ critical value of the parameter a (or d) is imaginary, 
so that there is no physical reversal speed and the wing under considera- 
tion is safe against reversal (in the speed range under consideration). 

If the critical value of the parameter a has the sign opposite to 
that of the value of ei^. (or the other value of eq used as a 

reference) or if the critical value of d has the sign opposite to 
that of the sweep angle A, the reversal dynamic pressure will be 
negative. In that case also the wing is safe against reversal, since 
a negative reversal dynamic pressure cannot be obtained at any real 
speed. 


Calculation of Control Power and Maneuverability 
The calculation of the twist distribution for a given aileron 


deflection is carried out in table Ill(b). The matrix 


is entered at the left and the column 




[l] - a [A]] 


at the right. Usually 



12 


MCA EM No. L8H24a 


it will te convenient to let 0^6 = 1 (except where modified hy 
fig. 2) and then miiltiply the resulting twist distribution and 

rolling power hy the true ccg6 values if desired. The colimm 

is then premultiplied hy the |^aJ matrix ^step or and 

entered in the second column at the right, which in turn is multiplied 
hy —a to yield the third column. The simultaneous equations with 
the coefficients at the left and the knowns at the right (the third 
column) are then solved for the unknown values. 

It will he noted that if the same values of a are selected as 
were used/ in the calculation of the aerodynamic loading hy the method 

of reference 6, the [[I] - a[A]] matrix will already he available. If, 

in addition. Grout’ s method of solving simxiltaneous equations has been 
used to solve the simultaneous equations, the auxiliary matrix will 
also he available, so that calculation of the ag values for the 

aileron loading will require very little time. 

In some of the special cases discussed in the preceding sections 
care must he taken to use the proper parameters in conjunction with the 
matrices calculated for these special cases. In the case where e^^ 

is zero, the a values must he based on the reference value of e^^ 
selected in calculating the matrix; in the case where e^^ is zero along 

the entire span the parameter d must he used instead of a in 
table- Ill(h) . 


The res^llting 
multiplied hy 


values may he added to the effective ctg6 
and plotted over the span to yield the 


values. 


/cc^X 

aerodynamic load distribution / ^ — j . The rolling-moment coeff i— 

V°“eV=i 

cient due to this forcing loading (over both wings) may he obtained 
from the relation 



This coefficient, which is a direct measure of the rolling power, is 
seen to be dependent only on q/qp (except for the factor m^), since 

— and an is consteint for a given speed range, 
qp ap 
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The rolling maneuverability depends not only on the rolling power 
but also on the damping in roll. The rate of roll per unit aileron 
deflection is given by 


/pb\ 

W6=l 


= 06 




(16) 


where . forcing coefficient ceilculated from equation (I5) 

with 005S = 1 and damping coefficient calculated from 

equation (40) of reference 6 for a value of o,g^= = 1 at the wing 


tip. 


Illustrative Etcample 

The method described in the preceding sections has been used to 
analyze the lateral maneuverability of the wing considered in the 
illustrative example of reference 6. The required additional parameters 
of this wing are presented in table IV (a), which follows the fdrm of 
table I. For convenience a value of o6& = 1 has been selected. The 

equivalent value of ct5& at the station ~ = O.h is obtained from' 

®w 

figure 2 for the given values of yi/s^ and yg/s^. The auxiliary 
elastic matrix for the subsonic case has been calculated by following 
the form of table 1 5 the re stilting matrix is shown in table IV (a). 

The aileron-reversal matrix for the subsonic case is calculated 
by means of the form of table II. Sevei^ of the steps, as well as 
the result, are shown in table IV (b) for the subsonic case. Iteration 
of the aileron-reversal matrix (by means of the form of table III (a) 
or otherwise) yields a value of aj^ = 2.36b. A similar calculation 
for supersonic speeds yields a value of a^ = O.I280. From these two 

values and the definition of the parameter a (see reference 6) the 
dynamic pressure required for reversal may be calculated and plotted 
against Mach number, as shown in figure 3. Also shown in figure 3 
for conqiaidson are the dynamic pressures required for divergence as 
well as the actual dynamic pressures at several altitudes. Where 
the dynamic pressure required for reversal is less than the actual 
dynamic pressure the aileron control is reversed. For the example 
wing reversal occiars at a Mach number of I.3, approximately, at sea 
level. 


( 
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The aerodynamic 
calculated "by means 


loading due to aileron deflection has teen 
of the form of tatle 111(1)). For the subsonic 


case and a value of 

shown in table X(b) 
of table Ill(b) are 
values of the final 


a = 0.552 the [M - a m] matrix is that 

of reference 6. The three columns at the rig^t 
given in table V. Also shown in table V are the 
matrix calculated from the third column, the 


values of 



the values of 



for the twist distri- 


bution, which are obtained by multiplying the values of the preceding 


two columns by each' other, and the values of 




for the 


aileron-deflection distribution, which are obtained by multiplying the 


058 


values by the 



values. 


The aileron-deflected distribution given by the last column applies 
directly only to the rigid- wing casej it is plotted as such in figure 4. 
It will be noted that the equivalent a58 value of figure 2 affords 
a convenient guide for fairing or rounding off the o§8 distribution. 
For the flexible wing the ceilculated twist distributions, such as the 
one shown in the next to the last column of table V, must be added 
algebraically to the aileron-deflected distribution. This is best 
done by first plotting them separately and then adding them point for 
point to the aileron distribution. The net distributions obtained in 
this manner for several cases are shown in figure 4. It is seen 

immediately that the distribution for case 5 (supersonic speeds, 

G[ \ \ 

— = -l.OOj indicates that the wing is operating at a speed above its 
reversal speed j actually the ratio for this case is I.I54. The 


moments of the twist and aileron- deflected distributions are obtained 
by multiplying them by the first row of the K2 niatrix. The rolling- 

moment coefficient is obtained from equation (I5) or by adding the 
moments of the aileron-distribution curve eind the twist curve 


algebraically and multiplying the result by 



The ratio 


of the flexible-wing rolling-moment coefficient obtained in this 
manner to the corresponding rigid- wing rolling-moment coefficient is 
plotted in figure 5(u) against the ratio q/qj). The lateral maneuver- 


ability is calculated by means of equation (I6) using the damping 
coefficients calculated in reference 6 and is also plotted as a fraction 
of the rigid- wing value in figure 5(a). It will be noted that both the 
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maneuveratility and the control power hecome zero at a value of 
~ = -0.87, which is indeed the ratio of the reversal to the divergence 
dynamic pressure at supersonic speeds, as is shown in figure 5. 

Since the ratio has "been determined as a function of altitude 

and Mach number in figure 3 the parameters of figure 5(a) can he 
plotted as functions of altitude and Mach number, as has been done in 
figure 5(1)" It is seen that the maneuverability and, to a lesser extent, 
the control power are relatively low at supersonic speeds, particularly 
at low altitudes. Since at high speeds even a small value of pb/2V 
implies a fairly large value of the rate of roll p, this sitviation is 
not necessarily alarming. The wing in question should have adequate 
control at all speeds for altitudes greater than about 20,000 feet. 


DISCUSSION 


The discussion of the aerodynamic and structural assumptions of 
reference 6 is pertinent to the analysis of this paper as well. The 
additional aerodynamic assumption made in this paper, to the effect 
that aerodynamic induction effects may be estimated -by reducing the 
strip- theory lift distribution by an over-all correction and rounding 
off the distribution both at the wing tip and at the aileron ends 
(using the equivalent values of fig. 2 as a guide) , is consistent 
with the other aerodynamic assumptions. The reduction of the load 
distribution appears to be the same as that for a linear twist. A 
more refined way of taking the induction effects into account would be 
to use aerodynamic influence- coefficient matrices. As soon as 
suitable aerodynamic matrices become available they may be included 
in the method of this paper. 

G?wo additional structural assumptions are made as well. In the 
first place, it is assumed that the angle 6 between the wing and the 
aileron is constant along the span. This assumption appears to have 
been made in almost all of the published Investigations into the 
problem of lateral-control reversal and appears to have yielded 
satisfactory results. The, shorter the aileron and the greater the 
number of points at which the aileron is supported and at which its 
hinge moment is taken out the more nearly true the assumption would be. 

In the second place, it is assiamed that the control linkage is stiff, 
so that the aileron angle for a given stick displacement is independent 
of the dynamic pressure. This assumption need not be nwriFt if it is 
kept in mind that the results . obtained by the method of this paper are 
for a given aileron angle and that the true aileron angle may be less 
at high dynamic pressures than at low ones. Thus, in order to account 
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for the control-linkage deflection, it is necessary only to calculate 
the ratio of the true aileron angle at a given dynamic pressure to that 
at zero dynamic pressure for the same stick position. The calculated 
control moment and maneuverahility must then he reduced hy this factor. 
Since deformations of the control linkage only affect the aileron 
effectiveness, they have no bearing on the reversal speed. These 
deformations may lead to aileron divergence for wings with heavily 
overbalanced ailerons. This problem, as well as the problem of wing- 
aileron divergence, has not been considered in the present analysis, 
however . 

The fuselage and tail do not contribute any appreciable amounts 
to either the control or the daarping moment, so that their effects may 
ordinarily be neglected for the puii)ose of lateral- control calculations. 
Similarly, the effect of wing camber does not enter into the problem 
because the only important effect of camber is to give the flexible 
wing a symmetrical lift distribution if it is set at the angle of 
attack which would give zero lift for the rigid wing; this symmetrical 
lift distribution has no effect on the lateral- control problem. 

As in reference 6, the effects of the inertia loading on the 
aerodynamic loading have not been considered explicitly in the analysis 
of this paper. As pointed out in reference 6, however, the structural 
deformations due to the Inertia loading may be calculated conveniently 
by means of the integrating matrices and then considered as part of 
the geometrical angles of attack. This procedure may be applied in 
the case of a rolling wing to determine the change in rolling moment 
for a unit rolling acceleration at any given Mach number and dynamic 
pressure. This rolling moment must be taken into account in estimating 
the rolling accelerations due to a given forcing moment at any time 
before the steady-roll condition is reached. 

At transonic speeds there is considerable uncertainty in the 
aerodynamic parameters. The control power is directly proportional 
to the value of the parameter cjg = which may be quite low 

in the transonic region due to the fact that the aileron is located 

in a region of separated flow. The method of this paper is applicable 
to this case if the value of cjg is known for the rigid wing. If 

the decrease in this parameter due to flow separation is 40 percent 
at a given Mach number and the loss in control power due to wing 

flexibility amounts to 20 percent, for example, then the total loss 

is 52 percent (l - O.60 x O.80). The loss in maneuverability due to 
the decrease in c^^ may be much less than the loss in control power, 

however, since a decrease in cjg would usually be accompanied by a 

decrease in m^, and hence in the damping coefficient. 


t 



NACA EM No. L8H2lfa 


17 


Should the value of the parameter decrease to zero or reverse, 

aileron reversal will occur. This type of reverseil is altogether 
different from the type of reversal discussed in this paper, since it 
is due entirely to aerodynamic action, whereas the reversal of concern 
in this paper is due to aeroelastic action. Both types of reversal 
are largely independent of each otherj aerodynamic reversal will occur 
at a given speed regardless of the stiffness of the wing, whereas 
aeroelastic reversal will occur ordinarily at a different speed which 
is unaffected hy the aerodynamic effectiveness. 

As pointed out in reference 6, the method on idiich the analysis of 
this paper is Based does not require the semirigid representation 
employed in many analyses of aerodynamic loading of flexible wings and 
of lateral- control reversal. The method of this paper takes the actual 
stiffness distribution and plan form into accountj it integrates the 
differential equations exactly (within the accuracy of integrating 
matrices) without simplifying the wing to one of constant- chord segments 
with all the flexibility concentrated at the ends of the segments and 
without requiring time-consuming graphical integrations. Exurthermore, 
the method of this paper furnishes the aerodynamic loading and hence 
the control power and maneuverability directly without iteration. An 
iteration is required to calculate the reversal speed, but this 
iteration is straightforward in application and converges rapidly in 
most practical cases. If it is preferred, the iteration may be 
dispensed with and the critical value of the parameter a determined 

instead by setting the determinant of the matrix [[I] - a[A]] 

equal to zero. This, procedure implies calculating the coefficients of 
and solving a sixth-degree or tenth-degree equation in a (depending 
on whether the 6-point or the 10 -point method is used), so that it 
is ordinarily more laborious than iterating the matrix. 

Some general effects of sweep and of the moment arms e^ and 62 

on the aeroelastic reversal speed may be of interest. The ratio of 
the reversal parameter aj^ of a given wing to that of the unswept wing 

obtained by rotating the given wing aj^ 'is shown in figure 6(a) 

plotted against a function of th© sweep angle for subsonic and supersonic 
speeds} the two curves were obtained by considering the wing of the 
illustrative example to be rotated in such a manner as to keep the 
®lj.Cr cos A (gj) 

parameters , as well as the chord, stiffness, and 

moment arm (e^ and distributions constant. ’ 

It appears that both sweepback and sweepiforward tend to decrease 
the reversal parameter and hence the reversal speed. At supersonic 
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speeds, or more specifically, at s m a l l values of the parameter 
eg^Cr cos A rr^ 

1/7—^ the reversal speed for the sweptforward wing is 

®A I (GJ Jr 

somewhat lower tlian that of the swepfbach wing, whereas at hi^er 
values of the parameter the variation of the reversal speed with the 


sweep parameter tan Ai 


(GJ)i 


is more nearly symmetrical with respect 


n (El)r 

to the zero- sweep case. There are some Indications that this behavior 
is not typical, of all. wings hut rather is due to the fairly large 
variatlcais of the e^ values, which also appear to he responsible for 

the deviation from linearity of the d - a c\irve of figure 7 of 
reference 6, as well as of the e2 and values over the span of the 

exan5)le wing. In general it appe ars tha t, for small values of the 

^ [JEIi 

^A \ (GJ)r'' 


moment-arm parameter 


®2^°r 


the variation of the reversal 


speed with the sweep parameter tan 


(C 3 J), 


should he nearly symmetrical 


, (El)j. 

and that, for large values of the moment-arm parameter, the reversal 
speed should tend to he lower for swepthack than for sweptfoivard wings. 


The variation of the reversal speed of an unswept wing with the 
moment- arm ratio is shown in figure 6(h) for wings which 1 have the 
same distributions of the parameter ©ny^e^^ and ©^©2^. ©long the span 

hut have different values of ei^ and 02^- The parameter aj^^ is 

plotted against the ratio — i — , where the value of € is selected 

1 + 

at the mid-aileron station. It is seen that the plot 1s linear for ho-th 
the subsonic and the supersonic case. The diff©rence in these cases 
is due to the different variation of 02^ and &2 along the spanj if 

the variations were the same or if e^ and oq were constant a1 nng 
the span the two lines of figure 6(h) would coincide. Since the 

reversal parameter ap is proportional to — i — , and since the 

^ 1 + 

reversal dynamic pres8^lre is directly proportional to the reversal 

parameter a and inversely proportional to the valiie ei (by 

■‘•r 

definition of the parameter a) , it follows that the reversal 

dynamic pressure is approximately proportional to the ratio — — ^ . 

+• ©2 

From figure 1 it is seen that the sum of e^ and eg represents 
the distance from the aerodynamic center to the center of pressvire of 



NACA EM No. L8H2l^a 


19 


the lift due to aileron deflection and is independent of the location 
of the elastic axis. This fact corroborates the commonly made 
observation (see reference l) that the reversal speed is independent 
of the location 'of the elastic axis in the case of unswept wings. 

The control power and maneuverability cannot be related to the 
structural geometric parameters in as relatively simple a manner as the 
reversal speed. The control power is a function of both the ratio 

and the ratio it normally decreases with the rate of 

decrease being slow at first and then more rapid for positive values 
of (which would generally be obtained for unswept and swept forward 

wings) and vice versa for negative values of (which would 

generally be obtained for sweptback wings) . The variation of the 
maneuverability should generally be similar to that of the control 
power, since the dan5>ing coefficient decreases (or in the case of 
unswept and sweptforward wings increases) steadily with q/qj) and is 

independent of ^ . 

From the calculations for the example wings it appears that the 
control power and maneuverability of sweptback wings tends to be relatively 
low, particularly at supersonic speeds. If it should happen that a 
combination of high sweep and large moment arm b 2 leads to an \mdeslrably 

low maneuverability and neither of these parameters can be changed, it 
may be necessary to resort to unconventional control devices. Leading-edge 
ailerons, for instance, have negative values of the moment arm B 2 , so 

that wings equipped with them tend to reverse at very high speeds, if at 
all. This configuration has the additional advantage of relatively high 
effectiveness at transonic speeds. On the other hand, the effectiveness 
of leading-edge ailerons at subsonic speeds is so low that they wooild 
have to be used in conjunction with trailing-edge ailerons if a great 
deal of flying were to be done at subsonic speeds. Furthermore, they pose 
a number of difficult structural and other design problems, so that it 
would be well to consider them only as a last resort. 

Another means of raising the reversal speed and of increasing the 
control power is to change the stiffness of the structvire. In general, 
lateral control can be improved by increasing the torsional stiffness 
or the bending stiffness. In some cases, however, the increase of the 

reversal parameter due to a change in the parameter tan A 

(see fig. 6(a)) produced by a decrease in the torsional stiffness (GJ)^ 

may be so rapid as to cause a net Increase in the reversal speed. Finally, 
any increase in the purely aerodynamic effectiveness of the aileron- 

airfoil combination results in a proportional increase in the lateral 
control effectiveness. 
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CONCLUDING REMAEKS 


A method has "been presented for calculating the effectiveness and 
reversal of lateral control as well as of the aerodynamic loading and 
rolling moment produced hy aileron deflection on fleiihle wings of 
arbitrary plan form and stiffness. 

It has been shown that the aileron-reversal speed decreases with 
both sweepback and sweepforward and that the effectiveness of conven- 
tional aileron configurations on sweptback wings at supei^onic speeds 
tends to be relatively low. The control effectiveness and the resulting 
maneuverability of the airplane may be increased by varying the structural 
stiffness and, if necessary, resorting to unconventional control devices, 
such as leading-edge ailerons. 


Langley Aeronautical Laboratory 

Natloneil Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE IV* CALCULATIOM FOR THE E)(AMPLE WINS AT SUBSONIC SPEED 
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Figure 4.- Loading of example wing with aileron deflected. 
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Figure 5.- Control power and effectiveness of example wing. 




